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Abstract 

An exchange-correlation energy functional E xc and the resultant exchange- 
correlation potential v xc (r) in density- functional theory are proposed us- 
ing orbital-dependent coupling-constant-averaged pair correlation functions, 
g rTa ' (r, r') for electronic structure calculations of atoms, molecules, and solids. 
These orbital-dependent g a<7 (r, r') fulfill the symmetric property, the Pauli 
principle and the sum rules. In the limit of uniform density g aa ' (r,r') are re- 
duced to the very accurate analogues of the electron liquid that are obtained 
from an interpolation between long- and short-range correlations involving the 
exchange corrections. The major contribution of f X c( r ) is given in the form of 
the Coulomb interaction with the exchange-Coulomb hole around an electron. 
The present theory not only guarantees local charge neutrality, but also repro- 
duces the exact asymptotic form of the exchange potential, v x (r) = —e 2 /r for 
finite systems. The present method of dealing with correlations, if properly 
applied to finite systems, can give even the asymptotic form of the correlation 
potential v c (r) of order r~ 4 as well as the van der Waals potential of order 
r -6 for large r. 
PACS numbers: 71.15.Mb, 71.10.Cal 
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I. INTRODUCTION 



The most important subject in density-functional theory (DFT) [1,2] is the exploitation 
of an accurate exchange-correlation energy functional E xc and the corresponding exchange- 
correlation potential i> xc ( r ) beyond the local-density approximation (LDA) and the general- 
ized gradient expansion approximations (GGA's). In a previous paper [3] we have proposed 
an orbital-dependent correlation energy functional E c [{y?i}{£i}] which is to be used by the 
optimized potential method (OPM) [4] in combination with the exact orbital-dependent ex- 
change energy functional E x [{<£>,}]. It comprises a direct and exchange pair of second-order 
perturbation-like terms constructed from Kohn-Sham orbitals and Kohn-Sham energies, but 
one of the two Coulomb interactions entering each term is replaced by an effective interaction 
v eff (r) which can in principle contain all the effects of third- and higher-order perturbation 
terms. From the accurate knowledge of long-, intermediate-, and short-range correlations of 
the electron liquid, we have defined such effective interaction v eff (r) for the electron liquid in 
order to substitute it for the orbital-dependent correlation energy functional above. These 
orbital-dependent exchange and correlation energy functionals depend implicitly on the elec- 
tron density n(r) through Kohn-Sham orbitals and Kohn-Sham energies. The application 
of such implicit functionals to DFT is founded on the basic assumption that the correlation 
potential v c (r) as well as the exchange potential v x (r) will be evaluated by solving an integral 
equation which relates i> xc (r) to 5E XC /Sifi(r) and 5E xc /5ei. The procedure for solving this 
integral equation has been called the optimized potential method (OPM) [4]. Actually, the 
OPM has been applied to the evaluation of the exchange potential v x (r) with great success 
[4,5,6,7]. 

Very recently, it has, however, been pointed out [8] that for finite systems the corre- 
lation potential v c (r) obtained from a direct application of the OPM to the lowest-order 
contribution of E c , i.e., the second-order perturbation terms constructed from Kohn-Sham 
orbitals and Kohn-Sham energies, in fact, becomes divergent in the asymptotic region for 
large r. The unphysical behavior above is caused by the presence of unoccupied orbitals 
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and unoccupied energies in the perturbation expression for E c . This implies that the OPM, 
though apparently most promising, might not be valid for the evaluation of the correlation 
potential v c (r). 

In this paper we propose an approximation method which enables one to avoid the 
above difficulty in the OPM in order to utilize the effective interaction of the electron 
liquid we have introduced for the evaluation of the correlation potential v c (t) in DFT. 
For this purpose we transform the orbital-dependent exchange and the orbital-dependent 
second-order perturbation-like correlation energy functionals given in the previous paper 
[3] into another expression written in terms of the orbital-dependent coupling-constant- 
averaged pair correlation functions, g(r, r') [9,10]. This transformation enables one to deal 
with E xc as an explicit functional of the electron density n(r) and the major contribution of 
v xc (r) is physically appealing since it is given in the form of the Coulomb interaction with 
the exchange-correlation hole around an electron at the position r. The remaining terms 
involving the functional derivative, 6g(r', v")/8n{r) probably will make a minor contribution 
to f xc ( r ) since the six-fold integrations over r' and r" may be expected to make only a 
secondary correction to the Hartree potential, giving a rather small and nearly uniform 
modification in the electron density n(r) of the system. Therefore we deal with these terms 
in the spirit of the local-density approximation (LDA), i.e., calculate these terms in the 
limit of uniform density and replace their density dependence by the local density n(r) of 
the system. It guarantees that the total i> xc (r) in the present theory is reduced to the LDA 
in the limit of uniform density. 

In section II we give an approximate expression for the orbital-dependent spin-parallel 
and spin-antiparallel coupling-constant-averaged pair correlation functions for use in the 
exchange and correlation functionals in DFT. We show that these pair correlation functions 
not only fulfill the inherent symmetric property as well as the requirement due to the Pauli 
principle, but also satisfy the basic sum rules. In section III we give an approximate form 
of exchange and the correlation potentials for practical use. The last section is devoted to 
conclusions and discussions. 
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II. ORBITAL-DEPENDENT COUPLING-CONSTANT- AVERAGED PAIR 

CORRELATION FUNCTIONS 



Let us start with the general expression for E xc (r) that Harris and Jones [9] have first in- 
troduced in DFT. This expression gives a clear physical meaning of exchange and correlation 



In Eq.(l), the coupling-constant-averaged pair correlation function g(r,r') is defined using 
the technique that the electron density n(r) in the real interacting system is maintained while 
the Coulomb interaction among electrons is adiabatically switched on as a perturbation. 
The function thus defined, g(r, r') exhibits general behaviors analogous to those of the 
usual pair correlation function g(r,r') as a function of |r — r'|. The difference of g(r, r') 
from unity is, however, generally reduced in magnitude compared with that of g(r, r'), 
which is noticeable particularly for small separations. This is because a lowering in the 
electron-electron interaction energy due to the presence of g(r, r') from its value in the 
Hartree approximation is partially cancelled by an increase in the kinetic energy which is 
caused by correlation and at the same time cannot be described in terms of T s , the kinetic 
energy of the reference non-interacting system with the same n(r) as the real interacting 
system, through the coupling-constant-averaged pair correlation function g(r, r') introduced 
in DFT. From Eq.(l) it is evident that exchange and correlation have the effect to reduce 
the contribution of the Hartree energy functional from short separations and the effect to 
enhance the contribution from long separations as a counteraction. Correlation between 
electrons not only gives the net reduction in the electron-electron interaction energy that is 
gained at the cost of an inevitable increase in the kinetic energy, but also plays an important 
role in the determination of the electron density n(r) of real many-electron systems that are 
for the most part stabilized under the influence of the presence of nuclei. Correlation in real 
systems gives rise to a change in the electron density n(r) in contrast with correlation in the 
uniform electron liquid. 



in DFT. 




(1) 
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In addition to the effective interaction v eff (r) we have also calculated with high accuracy 
the spin-parallel and spin-antiparallel coupling-constant-averaged pair correlation functions 
for the electron liquid, g a(7 (\r — r'|) and g a ' a (\r — r'|) . We have then ascertained that 
the difference of g aa (\r — r'|) from the Hartree-Fock correlation function gnpQr — r'|) is 
significantly reduced in comparison with that of g aa (\r — r'|) and that the difference of 
g a ' a (\r — r'|) from unity is appreciably reduced in comparison with that of g aa (\r — r'\). 
These reductions are due to the correlational increase in the kinetic energy, which is by 
definition involved in the coupling-constant-averaged pair correlation functions. It should, 
however, be noted that these coupling-constant-averaged pair correlation functions for the 
electron liquid are inappropriate for use in Eq.(l) since no orbital-dependent features can 
be involved in these functions. 

Consider the second-order perturbation-like expression we have proposed in the previous 
paper for the orbital-dependent correlation energy functional [3]: 

E = 1 y» U j^' | (ViVjHri2)\VaVb)(^b\Ve ff {ri2)\^j) _ {ipupj \v{rv2)\(p a ipb) {¥a¥b\Ve S {r 12) \ipjipi) 

i,j a,b I &i ~\~ £j &a &b &i ~\~ £j &a &b 

{<Pi<Pj\v{ri 2 )\<p a <Pb) = j j dr 1 dr2(p*(r 1 )(p*(r2)v(r 1 2)(pa(r 1 )(pb(r 2 ), (2) 

( ( Pi ( Pj\ v <ff( r n)\ ( Pa<Pb) = j j ^ri^r2^*(ri)^*( r 2)^i/'( r i2)^a(ri)v?6(r 2 ), 

where </?j, ip a and e a are Kohn-Sham orbitals and Kohn-Sham energies, respectively. 
In Eq.(2), v{r\2) = e 2 /|ri — r2 1 and the effective interaction v eff (r) is the real-space Fourier 
transform of the interaction v eff (q) we have defined by an exact expression for the correlation 
energy of the electron liquid as 

*< s*< ,w ,(l-/(p + q))(l-/(p'-q)) 



^ = 2^) E E /(p)/(p>(q) : 



q p,o-,p',ct' £ p £ p+q + £ p' £ p'-q 

x ivM - $**'V eff (-p + p' - q)) , (3) 

where /(p) is the Fermi distribution function at zero temperature and e p = h 2 p 2 /2m; f2 
being the volume of the system. In Eq.(3), long-, intermediate-, and short-range correla- 
tions arising from all the higher-order perturbation terms beyond the second are taken into 



account in the form of v eff (q). We have evaluated the effective interaction v eff (q) by making 
a sophisticated interpolation between long-range correlation in the random-phase approxi- 
mation (RPA) [11] and short-range correlation in the particle-particle ladder approximation 
[12,13,14] such that the corresponding exchange interaction v eff (— p + p' — q) and its feed- 
back effect on the direct interaction v eff (q) are allowed for in a self-consistent way. The new 
expression for E c given by Eq.(3) not only reproduces available most accurate numerical 
values [15] of the correlation energy within an accuracy of 0.5mRy. per electron throughout 
the entire region of metallic densities, but also has the merit of giving spin-antiparallel and 
spin-parallel contributions of the correlation energy separately. It is important to realize 
that the spin-antiparallel contribution occupies about 70% of the total correlation energy 
and the spin-parallel contribution does occupy as much as about 30% throughout the entire 
region of metallic densities. Probably, this ratio will apply to all the valence electrons that 
take part in cohesion of real metals. 

It is important to notice that Eq.(2) can be transformed into the correlation part of 
Eq.(l). For this purpose we shall start with the spin-dependent version of Eq.(l), which 
enables one to deal with the exchange and the correlation energy functionals, E x and E c , 
separately. 

= 5E //drdr- e \<yV (r,r') - D, (4) 

where the spin-parallel coupling-constant-averaged pair correlation function g aa (r, r') may 
be split into the Hartree-Fock and the correlational contributions as 

?°M=g&(Trf+%°(T,i>). (5) 

Then the exchange energy functional E x can be written in terms of the Hartree-Fock spin- 
parallel pair correlation function ggp(r, r') as 

where 
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^ F (r,r') = l- ' 



E^( r )^( r ') 



(7) 



n a (r)n a (r') 

This function is constructed from occupied Kohn-Sham orbitals alone. On the other hand, 
the correlation energy functional E c consists of the spin-antiparallel and spin-parallel con- 
tributions as 

E C = J2 (Er + E?) , (8) 



Er = \SS ^ e2 ^- (r/) (r-(r, r') - 1), (9) 



£r = i//dr^ ^y gr(r,0, (10) 



c 



where g a ' a (r, r') is the spin-antiparallel coupling-constant-averaged pair correlation function 
and g^ CT (r, r') denotes the correlational part of spin-parallel coupling-constant-averaged pair 
correlation function. An inspection of these two contributions of E c manifests that develop- 
ment of the Coulomb hole and further lowering of the Fermi hole at short separations both 
make a contribution to the correlation energy. 

Now we can easily transform the second-order perturbation-like expression given by 
Eq.(2) into the above spin-dependent expressions written in terms of g°" a (r, r') and g^°"(r, r') 
by paying special attention to a single Coulomb interaction entering each of direct and ex- 
change terms in it. The two correlation functions, g a ' a (r, r') and g^i^, r') may be identified, 
respectively, as 

r - a{r T , ) _ 1= 1 g"g c ' ^(r)^_ g (r / )y og (r)^- g (rQ 

^^■(rjTl-o-^r ) ^ Ei a -\- £j-<j £ao ^b-cr 

x // *i* 2 ^ (T ( r i)^- -(r2)^e#(ri2)^ -(ri)^- C7 (r 2 ), (11) 

^o-(r)77- (T (r / ) ^ . a ^ + £j CT £ acr E^tj 

x y y > ^ri^r 2 v?L( r i)^L( r 2)M r i2){^(ri)v?i CT (r 2 ) -^(ri)^(r 2 )}, (12) 



where the electron density with spin a is given by 



occ. 



Mr) = E IvMr) 



(13) 



Note that these correlation functions are constructed from both unoccupied and occupied 
Kohn-Sham orbitals and Kohn-Sham energies. It is important to notice that these correla- 
tion functions in the limit of uniform density, i.e., if Kohn-Sham orbitals and Kohn-Sham 
energies in Eqs. (11) and (12) are replaced by plane waves and free electron energies, re- 
spectively, they are reduced to the accurate analogues of the electron liquid because of the 
very definition of the effective interaction v eff (r) that we have introduced. 

It is evident from Eqs. (11) and (12) that the present theory satisfies the symmetric 
property inherent in the pair correlation function as well as the requirement due to the 
Pauli principle. 



It is important to realize that the well-known sum rule concerning the exchange- Coulomb 
hole is in fact exhausted by the Hartree-Fock exchange hole alone. 



This is evident from the orthogonality between occupied Kohn-Sham orbitals. As a con- 
sequence, the Coulomb hole due to correlation between electrons with opposite spin and 
the correlational change in the Fermi hole separately have to integrate to zero. In other 
words, correlation due to the Coulomb interaction occurs among electrons, for either of two 
spin orientations, such that local charge neutrality is maintained at every position r of the 
system. The present theory has the striking merit that the requirements above are fulfilled, 
as is evident from the orthogonality between occupied and unoccupied Kohn-Sham orbitals 
appearing in Eqs. (11) and (12). 



^'(r,r')=r'V,r), 



(14) 



£"(r,r)=0. 



(15) 




(16) 
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Jdv'n. a {v') (^(r,r')-l)=0, 
|rfrV(r')^r(r,r / ) = 0. 



(17) 
(18) 



III. EXCHANGE AND CORRELATION POTENTIALS 

In order to obtain the exchange potential v xa (r) we evaluate the functional derivative of 
E x given by Eq.(6) with respect to n a (r). The major contribution of v xa (r) is given as 

/*'^te(r,r')-l). (19) 
J |r — r I 

This is nothing but the Coulomb interaction with the bare Fermi hole around an electron 
located at the position r. In the limit of uniform density it is reduced to 3/2/z X(7 [16] where 
/i xcr denotes the exchange contribution to the chemical potential of the uniform electron 
liquid. The remaining term of v xa (r) is written as 

The functional derivative entering Eq.(20) cannot possibly be evaluated for arbitrary elec- 
tronic systems, but Eq.(20) can be analytically evaluated in the limit of uniform density. It 
amounts to be — l/2/i X(T [16]. It is not likely that Eq.(20) may produce such an important 
contribution to v xa (r) as is comparable to Eq.(19) because of the six-fold integration over r' 
and r". Then we approximate this term by — l/2/i xo -(n CT (r)) following the spirit of the LDA. 

A similar approximation method may be applied to the correlation potential f C(7 (r). The 
potential v ca (r) in the present theory consists of the major functional derivative of E c given 
by Eqs.(8), (9) and (10) with respect to n a {v) and the minor correction as 

«~(r) = / tf^=^(?-°(T, r') - 1) + J dr' e ^lg^(v, r') + A^K(r)). (21) 

The first term is given in the form of the Coulomb interaction with the Coulomb hole around 
an electron located at the position r and the second term in the form of the Coulomb inter- 
action with a correlation-induced change in the Fermi hole around the electron. From these 
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expressions it is then evident that the first and second terms have the effect to significantly 
reduce the contribution of the Hartree potential from short separations, which is important 
to the accurate evaluation of the electron density n(r) of the system. This is because the 
screening of each attractive nuclear potential by valence electrons is expected to be signifi- 
cantly reduced at short distances by the first and second terms in Eq.(21) and the electron 
density of n(r) may be somewhat enhanced in the immediate vicinity of each nucleus. On 
the other hand, it is not likely that the remaining terms of v ca (r) may make such an im- 
portant contribution as is comparable to the first and second terms because of the six-fold 
integrations involved. Then we approximate the remaining terms of v ccr (r) in the spirit of 
the LDA. The last term A/i co -(n CT (r)) in Eq.(21) denotes such an approximation to those 
terms which involve the functional derivatives, 8g a ~ a (r',r")/8n a (r) and 5g™(r', r")/8n a (r). 
The coefficient A should be determined such that in the limit of uniform dfensity the total 
correlation potential v ca (r) is reduced to the correlation contribution to the chemical po- 
tential of the electron liquid, /i ca . We think that the contribution of v ca (r) involving the 
six-fold integrations will form a background-like potential against the Hartree potential and 
plays a secondary role in the occurrence of a correlation-induced change in the electron den- 
sity n(r) of the system. Then we have treated that contribution following the spirit of the 
local-density approximation (LDA). 

For the practical application of the present scheme, it is a delicate problem how to choose 
the density parameter in the effective interaction v eff (r) we have borrowed from the electron 
liquid. The best choice of the density parameter will be such an optimization as gives the 
minimum value of the calculated ground state energy. Another choice will be the effective 
average density of valence electrons that take part in binding of molecules or in cohesion of 
solids; for metals it will correspond to the nearly uniform density realized in the interstitial 
region outside muffin-tin spheres. 

Now all quantities requisite for DFT are available. The exchange energy functional 
E x is exactly given in terms of the Hartree-Fock spin-parallel pair correlation function 
ggp(r, r') and the correlation energy functional E c in terms of the spin-antiparallel coupling- 
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constant-averaged pair correlation function g a ' a (r, r') and the correlation-induced spin- 
parallel coupling-constant-averaged pair correlation function g™(r,r'). Correspondingly, 
the exchange potential v xa (r) and the correlation potential v ca (r) are also given in terms of 
ggp(r,r'), g a ' (T (r,r'), and g° a (r,r'), though the minor contribution of f XC r(r) and v ca (r) is 
approximately treated in a similar way as in the LDA. We have thus avoided the difficulty 
that is encountered in dealing with the correlation potential v ca (r) according to the OPM 
by transforming the exchange and correlation energy functionals into the expression given 
by Eq.(4) and making some approximations. 

In the present scheme, the major contribution of t> X(T (r) and v ca (r) as well as E x and E c 
depends implicitly on Kohn-Sham orbitals and Kohn-Sham energies through the coupling- 
constant-averaged pair correlation functions g aa '(r,r'). The present v ca (r), when applied 
to finite systems, is not divergent for large r, however. This is evident from its expression 
(Eq.(21)) and the sum rules (Eqs. (17) and (18)). A much higher degree of self-consistency 
is required to perform the present scheme than in the LDA or GGA's since the coupling- 
constant-averaged pair correlation functions involved are orbital- and energy- dependent. 

IV. CONCLUSIONS AND DISCUSSIONS 

Between metals and non-metals there is an essential difference in long-range correlation, 
or equivalently screening at long distances, depending on whether the Fermi surface exists 
or not. We have borrowed the knowledge of long-range correlation from the electron liquid 
in which perfect screening at long distances is realized in the presence of the Fermi surface. 
In this sense, the effective potential v eff (r) in the previous paper [3] is properly defined for 
metals. If one performs a similar interpolation between long-range correlation in the RPA 
and short-range correlation in the particle-particle ladder approximation for finite systems 
or semiconductors and insulators in order to properly define the effective interaction v eff (r) 
in the second-order perturbation-like correlation energy functional, the resulting effective 
potential v eff (r) takes the form of e 2 /e r at large distances where e is a constant of order 
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of the static dielectric constant. On the other hand, there is no essential difference in 
short-range correlation, or the analytic form of v eff (r) at short distances between metals 
and non-metals. The particle-particle ladder approximation [12,13,14] gives the adequate 
description of short-range correlation as well as of the cusp condition [17,18,19] on the many- 
body wavefunction, independent of whether the Fermi surface exists or not. Therefore the 
present interpolation method of dealing with correlations can give the correct asymptotic 
form of the correlation potential v ca (r) of order r -4 for finite systems as well as the van 
der Waals potential of order r" 6 [20], if the effective interaction v eff (r) for finite systems is 
properly defined as e 2 je^r for large r. 

From experiences it has been widely observed [10] that the LDA has the tendency to 
underestimate the nucleus-nucleus separation in the formation of solids and the tendency 
to overestimate the binding energy of solids. Probably, these shortcomings of the LDA 
will be, for the most part, ascribed to that the screening of each nucleus or ion by valence 
electrons is overestimated at short distances and that the ground state energy of the atomic 
state is evaluated too high in comparison with that of the solid state. In our view, the 
adequate inclusion of short-range correlation between valence electrons may correct the 
underestimate in the nucleus-nucleus separation in solids. It is important to notice that 
short-range correlation between valence electrons has the effect to give less screening of each 
nucleus or ion at short distances and therefore may significantly enhance the electron density 
n(r) in the immediate vicinity of each nucleus. As a counteraction of this enhancement, the 
nucleus-nucleus separation is expected to be somewhat increased. 

The present approximation to E xc in its construction is expected to provide a much 
more accurate evaluation of the ground state energy of both atoms and molecules or solids 
than the LDA and GGA's. In the present E xc , the self-interaction errors arising from the 
Hartree energy functional are completely cancelled by the Hartree-Fock spin-parallel pair 
correlation function ggp(r, r'), though the corresponding exchange potential v xa (r) involves 
a minor approximate part. Furthermore, the correlation energy functional E c in the present 
scheme cannot be overestimated in magnitude since we borrow the knowledge of the effective 
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interaction v eff (r) from the electron liquid instead of the density dependence of the correlation 
energy. The LDA overestimates the correlation energy between nearly localized or tightly- 
binding electrons due to the occurrence of the locally high-density region caused by these 
electrons since the correlation energy of the uniform electron liquid increases its magnitude 
like lnr s in the high-density region where r s is the usual density parameter. 

A clear understanding of the difficulty [8] that is encountered in the evaluation of v ccr {r) 
according to the OPM will require a detailed investigation of how to make OPM integral 
equations compatible with Kohn-Sham equations, i.e., how to combine the two different 
self-consistencies in a unified way. Probably, the unphysical behavior of v ca (r) reported will 
indicate that the authors may have attempted to combine the two self-consistencies in a 
forced way. It will be closely related to the fact that the integral equation for v xca (r), if 
treated in a way consistent with Kohn-Sham equations, involves the functional derivative 
5v xca (r)/5n a >(r') in it. The OPM integral equation in its exact form is equivalent to the 
Sham-Schluter equation [21]. This equation requires that the Kohn-Sham Green's func- 
tion G KS (r,r',e) defined for the Kohn-Sham reference non-interacting Hamiltonian and the 
standard many-body theoretical one-electron Green's function G(r, r', e) for the same many- 
electron system give the same electron density n(r). This requirement can be fulfilled if and 
only if both G KS (r, r', e) and G(r,r',e) are exact. According to a systematic perturbation 
expansion of the many-body Hamiltonian in which the Kohn-Sham noninteracting reference 
Hamiltonian is dealt with as the unperturbed system, the n-th order perturbation energy 
for n > 2 has turned out to be a functional of not only occupied and unoccupied 
Kohn-Sham orbitals and Kohn-Sham energies, but also of v xca (r) itself [4]. Therefore the 
n-th order integral equation for v xca (r) involves the functional derivative 5v xca (r)/Sn a '(r') 
and is extremely difficult to solve [4], if one goes beyond the first order perturbation E£) 
equivalent to E x . 

The correlation potential v ca (r) is by definition a functional of n CT /(r'). It is very proba- 
ble that the accurate spatial dependence of v ca (r) cannot be obtained without the detailed 
knowledge of its functional derivative, 5v xca {r) / 5n a i{r'). The situation is analogous to the 
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fact that the p- and u- dependence of the self-energy operator £ CT (p, uS) essential to the 
accurate evaluation of the quasi-particle energy dispersion E(p) for the electron liquid can- 
not be obtained without the detailed knowledge of its functional derivative with respect 
to the one-electron Green's function Gv(p', u/), 5£ CT (p, u)/5G a >(p' , u/) [22,23,24]. This is 
because the variation of £<r(p, ou) with respect to p or uj occurs only through that of each 
constituent in it, Gv(p',u/) [25]. Similarly, the accurate r-dependence of v xca (r) will not be 
obtained without the detailed knowledge of 8v xca (r)/8n a >(r') = 5 2 E xc /5n a (r)5n a i(r'). Prob- 
ably, the situation in DFT will be much more complicated than in the standard many-body 
formalism. This is because the simplicity of Kohn-Sham equations is established in DFT 
at the cost of the extreme complexity of how to obtain E xc and t> XC(7 (r). From a compari- 
son with the second-order functional derivative of the Hartree energy functional, e 2 /|r — r'| 
it can be seen that the factor of e 2 (g aa '(r, r') — l)/|r — r'| entering the present E xc will 
form a physically appealing and important part of the second-order functional derivative, 
5 2 E xc /5n a (r)5n a ,(r'). 

Finally we would like to stress that the present theory gives an adequate description of 
both the orbital-dependence and the non-local dependence of the basic two quantities E xc 
and v xca (r) in DFT through the orbital-dependent g aa '(r, r') we have constructed with the 
help of the knowledge of long-, intermediate-, and short-range correlations of the electron 
liquid at metallic densities and that it is a feasible scheme for practical calculations. 
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